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Abstract 

We study the behaviour of lightcone wavefunctions, in coarse transverse lattice gauge theory, 
when one or more parton lightcone momenta are small. This probes the limit of hadron structure 
at large and small Bjorken x. Finite-energy boundary conditions on boundstates allow one 
to derive the analytic form of wavefunctions in this region. This leads to simple, universal 
predictions for the behaviour of quark generalised parton distributions near their endpoints. 
For the first few meson wavefunctions at large N c , we give the simplest ansatz that incorporates 
all the boundary conditions. 
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1. INTRODUCTION 



Transverse lattice gauge theory |l|, |2| is a promising approach for ab-initio calculations 
of the hadronic light-cone wavef unctions. All successful calculations so far have been per- 
formed on coarse transverse lattices with effective 'colour-dielectric' degrees of freedom, 
keeping the longitudinal coordinates x~ = (x° — re 3 ) /a/2 continuous. Although this has 
required some parameters in the effective Hamiltonian to be fit phenomenologically, it 
has the great computational advantages that hadrons are dominated by states contain- 
ing few partons and there is a clean separation between the longitudinal and transverse 
dynamics. In fact, in this case, it is not much harder to analyse the dependence of wave- 
functions on the longitudinal coordinates than in 1 + 1-dimensional gauge theories. While 
the full wavefunction of a boundstate typically requires — even in 1+1 dimensions — a nu- 
merical solution, the behaviour when parton lightcone momenta k + vanish can be largely 
understood analytically. The prototype example is 1 + 1-dimensional large- N c QCD [3(, 
where the behaviour of the meson wavefunction for vanishing quark k + was found exactly. 
A similar analysis can be performed to some extent even in 3 + 1-dimensional large- N c 
QCD 0]. The key strategy in all cases is to demand cancellation of singularities in the 
lightcone boundstate Schrodinger equation when parton lightcone momenta vanish. Typ- 
ically, wavefunctions do not simply vanish in this high free-energy limit, as they would at 
large transverse momentum. Rather, there occur 'ladder' relations between states with 
different numbers of partons that couple the analytic structure of the infinite tower of 
lightcone wavefunctions. 

In this paper, we analyse in a similar vein the corresponding lightcone boundstate 
Schrodinger equation for coarse transverse lattice gauge theory at large- N c . By demand- 
ing finiteness of the boundstate energy, we determine the form of meson wavefunctions at 
small lightcone momentum. These depend on two exponents, a and j3, whose values are 
independent of all details of the wavefunctions, which encode the leading non-analytic 
behaviour. In general, naive intuition based on the free lightcone kinetic energy is wrong. 
In particular, wavefunctions become more singular as more momenta vanish. Quark dis- 
tributions diverge as x — > and all Fock sectors contribute as x — *■ 1. 

In the next section, we briefly review the colour-dielectric formulation of transverse 
lattice gauge theory and set up the boundstate Schrodinger equation to lowest non-trivial 
order of the colour-dielectric expansion that is used in the following two sections. Sec- 
tions|21and|3]determine systematically the behaviour of lightcone wavefunctions of mesons 
at small parton lightcone momentum. At first, in sectional we use the approximation 
where up to one gauge link- variable is allowed in a state. Then, in section 0J we relax this 
approximation and generalise many of the results, giving a simple ansatz that incorpo- 
rates all the boundary conditions. Applications of the results to the behaviour of quark 
generalised parton distributions are given in Section |5] and conclusions are discussed in 
Section 
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2. REVIEW OF TRANSVERSE LATTICES 



We introduce continuum light-cone co-ordinates x ± = (x°±x 3 )/v^2 and discretize the 
transverse coordinates x = (a; 1 , a; 2 ) on a square lattice of spacing a. Lorentz indices /i, v 
are split into light-cone indices a, (3 G {+, — } and transverse indices r,s G {1,2}. Quark 
fields ty(x + , x~, x) in the fundamental representation and longitudinal gauge potentials 
A a (x + , x~ , x) in the adjoint representation of the gauge group SU(N C ) are associated to 
the sites of the transverse lattice. Flux link fields M r (x + , x~ , x), which we choose to be 
in a complex matrix representation of SU(N C ), are associated with the directed link from 
x to x + ar, where r is a unit vector in direction r. In the transverse continuum limit 
a — > 0, these would be related to transverse gauge potentials A r (x + , x~, x). Subsequent 
analysis is done to leadingorder of the 1/N C expansion, which allows one to drop the site 
argument x of the fields jll . 



A. Lightcone Hamiltonian 

For finite transverse lattice spacing a, the Lagrangian may contain any operator invari- 
ant under transverse lattice gauge symmetries, Poincare symmetries manifestly preserved 
by the lattice cut-off, and renormalisable by dimensional counting with respect to the 
continuum co-ordinates x a . One may control the subsequent proliferation of operators 
either by approaching the transverse continuum limit, in which case only renormalisable 
operators with respect to co-ordinates x survive, or bv'strong-coupling' methods on a 
coarse lattice, such as the colour-dielectric expansion [5J. All calculations to date have 
used the colour-dielectric form and it is this approach that we analyse here. This 
involves fixing to the light-cone gauge A_ = 0, classically eliminating non-dynamical 
fields A + , 7 _1 3/, then expanding the resulting lightcone Hamiltonian in powers of the 
remaining dynamical fields M r , 7 + \I/. Truncation of such an expansion is a valid approx- 
imation provided wavefunctions of interest (typically those of the lightest eigenstates) 
are dominated by few-body Fock states. This is achieved by the choice of sufficiently 
large masses for the dynamical fields. When the remaining couplings in the truncated 
hamiltonian are tuned to optimise symmetries broken by the lattice cutoff, the choice of 
large masses corresponds to choice of a coarse lattice spacing a. 

The Lagrangian density L we consider, which was previously investigated in Ref. 0], 
contains all allowed terms up to order M 4 and \1/M\I/. (There are additional terms 
involving multiple traces of products of M fields that contribute only to the glueball 
sector at large- N c (if.) 

L = ~2??2 Tl { Fa/3F ^} + Tr {[(^ + i40 M r - iM r A a ][h.c.}} 

a,fi=+- r=l,2 

-fit Tr{M r M}} - V[M] + i^ a (d a + L4 a )* - fi f ^ 

+m a (^7 r M r * - ^7 r M r f *) + k s (*M r * + , (1) 

where F a/3 (x) is the continuum field strength in the (x°,x 3 ) plane and the link-field 
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potential is 



V[M] = -^^Tr{M r M s MtM]} + ^i-^Tr{M r MtM r Mt} 



a?N c 
'a?N c 



r=£s 



a 2 N, 
A 



Tr { M r M r Ml Ml } + ^ Tr | M * M a M^ M a > j . (2) 



<T=±2,cr'=±l 



We have denned M r = Ml„ and hold G\/N C1 K a ^/N ci and K R \/N r finite as iV r — > oo. 



In the chiral representation, 



2 l/4 



/« + \ 

'<'+ 



(3) 



decomposes into two-dimensional complex spinors 











-1 




1 H 


i: 











(4) 



The helicity subscript h — ± denotes the sign of the eigenvalue of 7 5 . Using Pauli 
matrices, we may then write various bilinears 



1 

71 



v/2 
i 

71 



(v)a 3 v + nV 3 «) 



ftru — V/V 



) 



V2 

We also introduce the matrices ( x which act on the 2-spinors: 



(5) 

(6) 
(7) 

(8) 
(9) 
(10) 



A 


1 


2 


-1 -2 


C A 


a 2 


-a 1 


-a 2 a 1 



11] 



In light-cone gauge v4_ = 0, A + and v are non-dynamical (independent of light-cone 
time x + ) and their constraint equations of motion are used to eliminate them at the 
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classical level. The lightcone Hamiltonian may then be obtained from the action (JTJ in 
the standard way in terms of the remaining dynamical fields u, Mi, and M 2 . If we define 

F = -u + ^mJ^ + I^Au, (12) 

J+ = uw t + i M x d- M{ , (13) 

A 

the lightcone Hamiltonian one finds is 

Under certain reasonable assumptions 0], the Hamiltonian (|14J) is a truncation of the 
most general Hamiltonian to order M 4 and u 2 M. It also contains some, but not all, 
allowed operators at order u 2 M 2 and w 4 ; most importantly, it contains those responsible 
for confinement [l||. 



(14) 



B. Fock space 



Introducing the mode expansions 

1 r°° dk~ 



M r (x 



U(X 



0,x~ 
0,x' 



4tT JO 



a_ r (/c + ) e 



+ al(k + ) e 



(^£; + ) e" ifc+a; " + t7V(fc + ) e ik+3 



we have 



aA,i# + ),a*, fc/ (^ 
a A ,jj(/c + ), a pM (k^ 
b h ,(k + ),b* h ,/~k+ 



0. 

SijShh' - k + ) , 
. 



(15) 

(16) 
(17) 
(18) 
(19) 



where i,j are colour indices, A and p G {±1, ±2} denote the orientations of link variables 
in the (a; 1 ,^ 2 ) plane, and (a\ t ij)* = (o>\)ji- The sequence of orientations {A 2 ,---A n _i} 
of link variables together with the P + momentum fractions x a = k^/P + and quark 
helicities h, hi are sufficient to encode the structure of Fock states contributing to the 
boundstate at zero total transverse momentum P = 0. Only gauge singlet combinations 
under residual gauge transformations in = gauge can contribute to states of finite 
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energy 0, so a boundstate containing quarks at large N c decomposes in terms of Fock 
states as 

00 / p+\ n/2 r p+ ( 

i^( p+ )> = E E 2 ^ ¥c (w) / tei—^H 1 - 

n=2 /i,A 2 ,...,A„_i,/i' V c/ JO y 

x ^ n (xi, . . . , x n ; h, A 2 , . . . A n _i, ti) 

xb{( Xl P + ) a\ 2 {x 2 P + ) ■ ■ • 4_ i (x n _ 1 P + ) dl,{x n P + ) |0) . (20) 
The state is normalised by 

i = E E / ■ ■ ' dx ™d 1 1 ~ E Xa ) i^ n ' 2 ( 21 ) 

n=2 h,\2,-,\ n -l>h' V a=l / 

if 

&(P+)\i/;(P 2 + )) = 2P+(27r) *(P+ - P+) . (22) 

The object of this paper is to determine the behaviour of the wavefunctions i/) n (xi, . . .) 
when one or more x a are small. The Fock space matrix elements of the invariant mass 
operator 2P + P~ resulting from (J 14)) were given in Ref. 0]. For ease of reference and 
to correct some typing errors in interaction (ii), we reproduce these in a slightly more 
efficient notation as Fig. Eland Table ITT1 in Appendix 1X1 Projecting 2P + P~|^(P + )) = 
Ai 2 \tp(P + )) , where M. is the invariant mass eigenvalue, onto individual Fock states, one 
obtains integral equations for the wavefunctions. 



3. ONE-LINK APPROXIMATION 

To begin analysing the small x behaviour of the boundstate equations, we make the 
approximation of restricting the sum in Eq. (j2*U|) to n < 3, i.e. at most one flux link in a 
state. We correspondingly restrict to interactions (i)(ii)(iii)(vii)(viii) (see Appendix IA"|) . 
The 7r or rf meson state can then be written as 

|^(P + )) = 2,/^P + / dx^(x)b\xP + )a 2 d\(l-x)P + )\$) 
V N c Jo 

o r- ( p+\3/2 pi rl-x 

+ V \ T ' / dx d y y^U s {^y)bKxP + )a 2 a\{yP + )d\{l-x-y)P + ) 
+^ A (x, y) b\xP + ) C V a{(yP + ) d\(l -x- y)P + ) } |0> . (23) 



To obtain the ^ = spin projection component of the p meson, one replaces a 2 — > a 1 
in the above. The wavefunctions of the Jj, = ±1 components of the p meson are more 
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complicated: 

|^( J P + )> = 2. — P + / dxi){x)b\xP + )^d\{\ -x)P + ) |0> 

V N c Jo 

2v/?(P+)3/2 fl dx f X dyY,{M^y)bKxP + )a"a\{yP + )d\{l-x-y)P + ) 



o Jo 



y) b\xP + ) C V 4(yP + ) d\{\ -x- y)P + ) 
+^ A (x, y) b\xP + ) ^CVV 3 4(yP + ) dt((i - x - y )p+) ) |0) . 



(24) 



This is the state that is odd under reflections V 1 : x 1 — > —a; 1 ; the even state is obtained 
by removing the a 3 . 

The bound state equations for these four states are almost identical. The n = 2 bound 
state equation is 



M 2 ^{x) 



Hx) + ^ + ^ N < 



x(l — x) 
G 2 N C f 1 dy 



dy 



l-x 



dy 



2tt J (x-y) 2 



o X V Jo (! - x )y. 

(ij}(y) -ip(x)) 



ip(x) 



"K 



d ^ {2x-y)jj s {x-y,y) | f 1 x dy (2-2x-y)il> s (x,y) 



7T 



l-x 



y/yx{x-y) 
2 dy ^^ A< " X ~ V ' V ^ 



y/y{l - x){\ - x - y) _ 



x(x — y) 



' L dy (l-x)(l-x-y) 



i^a{ x i y) f° r pi° n 

-iPa( x , y) for rho, J 3 = 
-ip' A (x,y) for rho, J 3 = ±1 



(25) 



For each of the wavef unctions, x/ig, ipA, V'aj ^ ne n = 3 bound state equation is 



M 2 ij(x,y) = ptj ( - + - 



47T 



H i>(x,y) 

x-y J y 

x +y y + z 

dz — =; — ip(x + y — z,z) 



o 



l-x 



dz 



yz(y - z) 2 

y + z 

yz(y - z) 2 



ip(x, z) 



x+y dz r 1 -* dz 
+ 



-2ijj(x, y) - - . 

\y Jo {y-z) Jo [y-zr 



(26) 



For ips{ x ,y), we include 

HfK s \fW c 



nxy 



2x + y 
x(x + y) 



ip(x + y) + 



2-2x-y 



(1 - x)(l - x-y) 



ip(x) 



(27) 




FIG. 1: The n = 3 wavefunction ip{x, y) in the limits x — > 0, y — > 0, or p = x + y — ► 0. 



For iI)a{xi y) in the pion and the rho ^3 = state, we include 

HfKasfN^y tp(x + y) 1 I ij){x) for pion 

v^F x(x + y) (l-x)(l-x-y) \ for rho 

For ipA(x,y) in the rho J7~3 = ±1 states, we include 



(28) 



| lifK a y/N~yj){x + y) 
^/txx[x + y) 

For ip' A (x,y) in the rho = ±1 states, we include 

HfK a y/N c yif)(x) 



0F(1 -z-y) ' 

Next, we analyze the singularities as one or more of the momenta go to zero, 
consider each of the following limits: 

lim^(x) = C2X a 
lim tp(x,y) = x a g(y) 

x— >0, 
y constant 

lim ip(x,y) = y p g\x) 

y—fO, 
x constant 

lim ip(x,y) = (x + yff 

x/y constant 



X 



x + y 



p b f(w) 



where x = wp, and y = (l — w)p. See Figure^ The remaining limit of interest, 1 
with (1 — x)/y constant, is equivalent to the limit (JH3|) by charge conjugation. 



(29) 

(30) 
We will 

(31) 
(32) 

(33) 
(34) 
■x -> 
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A. n = 3 equation, limit x — > 0: 

For y), we obtain: 

„_1 HfKsVK 



= ^ gdy)*?' 1 - 

G 2 N r f x+y di 



47T 



y + z 



ijj{x + y — z) — 2ip(x, y) 



(35) 



For tl>A(x,y), we obtain: 



= n) g A {y) X- 1 + ^""J^S fo)*- 1 



g 2 ^ dz 



47T 



y + z 



ijj(x + y — z) — 2ip(x, y) 



(36) 



In order to cancel the divergence of second term, we must set a = 0. We then find that 
the integral term, which is also 0(x a_1 ), goes to zero. Simplifying, we get: 



K s \/W c 

limips(x,y) = J- ^{y) 



lim iPa{x,v) 

x— >0 



KaVN c 



(37) 
(38) 



B. n = 3 equation, limit y — > 0: 

For this limit, the resulting equation is the same for ip$ an d ipA because the differing 
terms are non-leading. 







i4 



G 2 N r 



71 



g\x) /- 1 



G 2 N, 



4:71 

+0(y 



c _x 

2/ 



V.'/ 



dq 



1 + g, 
v 7 ^ 



(z) - 2y*y (z 



(39) 



We find the constraint {3 < |, so that the integrals converge, in which case the 0(y - ^) 
terms are non-leading. We also have the constraint (3 > 0, so that the expectation value 
of 2P + P~ is finite. Thus, we obtain: 







G 2 N r 



IT 



fo) y^ 1 - ^ M y p ~ l 2 M tan W + 2 ] • ( 4 °) 
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Simplifying, we derive: 



0=W- G 2 iV c /5tan (tt/3)] g'(x) y^ 1 . 



(41) 



Note that g'(x) is not fixed by this analysis. Eq. (J41j) is the same expression that one 
obtains for the two-link bound state equation without quarks 0]. 



C. n = 3 equation, limit x + y — > 0, x/y constant 
In this limit, we obtain, for both fs{w) and /^(tu): 







1 — w 



G 2 N, 



Air 



dq 



o (q-wy 



(2-w-q) f(q) 



2/H 



47T 



o (9-^) 



y/(l-w)(l-q) 
(w + q — 2qw) f (q) ( w\ a 



y/(l-w)(l-q) \q 

for 



2w/(«;) 



+0 {p' 1 ' 2 ) . 



« a \/i - w for 



(42) 



In order for the last term in the equation to be cancelled in the limit p — > 0, we must set 
b = a — |. Thus, 

lim ^(ar,j/) =P a-1/2 /(w) • (43) 



x+j/— >o 



D. n = 2 equation, limit x — > 0: 

Taking the x — > limit of Eq. (f2~K|) and using p3|l. we obtain 





G 2 N C C 2 


+x a " 1 / 
Jo 


dz 


~k 2 s N c C 2 


l-z 


7T 


+x a ~ 1 [ 
Jo 


dz 


~K 2 a N c C 2 


l-z 


7T 


+0(x°) . 







„a-l 



k 2 s N c C 2 ^ f K s ^W c f(2-z) f s (l - z) 



7i 



7T 



'■{V~zf A {l-z)) 



(44) 



For the expectation value of the energy to be finite, we must have < a < 1. 
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w 

FIG. 2: /sH and f A (w) vs. tu for couplings n f = Q.m2G*JN~ c , fi b = 0.2G^fW c , and a = 0.4. 
E. Self-consistent solution 

To determine a, one must find fs(w)/C2 and /a(w)/C2 that self-consistently solve 
Equations (|4"2|) and (|44jl. We obtained solutions numerically by discretizing w on the 
interval [0,1]. Using numerical calculations of fs(w)/C2 and /a(iu)/C2, one can study 
the dependence of a on the couplings in the Hamiltonian; see Figures OHU Note that 
a is very nearly constant with respect to link field mass fib (Fig. Also, note that 
a is outside of its allowed range for some areas of the coupling space, corresponding 
generally to large values of hi s \fW~ c and K a \TN~ c (Fig. HJ. Thus, the endpoint behavior of 
the wavefunction can restrict the allowed values of the couplings. 

4. BEYOND ONE LINK 

When relaxing the truncation on n and including the other interactions of Fig. 
most of the complication occurs in the analytic auxiliary functions, such as /, g, g' in 
the previous section. If we wish to determine only the form of the leading non-analytic 
behaviour of wavefunctions at small momenta, the analysis is somewhat simpler. We note 
that this is independent of spin, represented by quark helicity and link field transverse 
orientation. Consider the effect of including n = 4 wavefunctions ^4 and the other 
interactions of Fig. El In the following we ignore the spin dependence, finding only the 
momentum dependent part of these contributions. Accordingly, we label wavefunction 
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FIG. 3: a vs. fj? f and k s = -0.323G, K a = 0.162G. 




2.5" 



FIG. 4: a vs. k s /G and K a /G; fi f = S.m2G^fW c , and fi b = 0.2G^ C . 



components ip n by the number of partons n and their lightcone momenta only. 



A. n = 2 equation 

To the n = 2 boundstate equation (|2~5j) we must add interactions (iv)(vi): 



dydz 



y-z 



4iT(y + zfJyz 



— V — -2? V, z ) + antiquarks 



from (iv), and 



dydz 



47r(x — z)y/yz 



ip4(x — y — z, y, z) + antiquarks 
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from (vi). If we assume that these contribute at the leading x a 1 order as x — > in the 
analsyis of Section 13 Dl then 

ij) A (x,y,z) ~ (1 -x) ' 1 as x -> 1 . (47) 



B. rt = 3 equation 

In the n = 3 equation (|2bj). we must include interaction (u), that connects n = 3 to 
n = 3 and was previously omitted: 

/ete — t r i>3(x + y — z, z) + antiquaries . (48) 
47r(x + y)y/yz 

This is subleading in the x — > analysis of Section 13 Al and the ?/ — > analysis of 
Section [3 Bl (recall that (3 < 1/2). It contributes at the leading order p b ~ x in the analysis 
of Section l3~C^ but does not affect the relation b = a — \. To the n = 3 equation (J2*HJ) we 
must also add terms involving ^4 via interaction (ii): 

/ dz — ^= (- ± — ?— ) ^ 4 (z - z, z, y) . (49) 

7 2^7TZ \X X — Z J 

If we assume that these contribute at the leading x~ x order of the analysis in Section 15 Al 
then 

ipi(wp, (1 — w)p, z) ~ p~ l l 2 as p — > . (50) 

In the y — > analysis of Section 13 Bl ()49|) does not contribute at leading order since, by 
normalisability, /3 — 1 < — 1/2 and ^4(2; — z, z, y) cannot diverge faster than 1/ y/y when 
y — > 0. In the analysis of Section PTO property (|4T|) implies that (J4T?|) contributes at the 
leading p a ~ 3 / 2 order as p — > 0. 



C. n = 4 equation 

Finally, there will be a new n = 4 boundstate equation for Ai 2 tp4 1 (x,y, z). The mo- 
mentum dependent parts of these contributions are given in Table |U 

Repeating the analysis of Sections l3 Al and l3 Bl for n = 4, we again find from comparing 
interactions (i) and (ii) that 

ip4(x, y, z) ~ x° as x — > , (51) 

since all other interactions are subleading, while from interactions (viii)(ix) 

fa(x,y,z) ~ y p as y -> . (52) 

We confirm directly, from the n = 4 amplitudes of Table I, the correctness of Eq. (I5UI) 
and Eq. (}4*Tj) . thus verifying the assumptions that went into deriving them. In a similar 
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TABLE I: The momentum dependent parts of interactions contributing to the boundstate 
equation for Ai 2 ^^, y, z). 



-ipi(x, y, z) + antiquarks 


(*) 


TTfTT ( r ± ) + U,z) + antiquarks 


(it) 


{V l$tlt- V m + antiquarks 


(iu) 


/ du ^(a+v) /w } + anti q uarks 


(«) 


t 2 (x+y)% + antiquarks 


(vi) 




(vii) 


^ p j du (^^w) + antiquarkg 


(viii) 


-lp (u+z)(2y+j;-«)'!/>4(a:,y+z-D,u) 
8,r J -y/yzi;(j/+.2-u)(.z-?;) 2 


(ix) 


yzv(y+z—v) 


(x) 


07r J yjyzv(y+z-v)(z+y) 2 


(xiii) 



way we also deduce as p — > 0: 



and 



^4(2;, (1 — ~ p° (53) 

^ 4 (wp, y, (1 - ~ p p (54) 

i/j 4 (wp,y,l-y-p) ~ p ° (55) 

^(z,y,z)~(l-y) a_1/2 aflj/->l (56) 



D. Summary 

We can summarize all these results in the following forms for ip n , which have the 
correct endpoint behaviour when any number of momenta vanish: 

fo(x) = C 2 x a (l-x) a 

yP(l - y\<* 

= C3 [(x + y)(l-xW-^ (5?) 



ip 4 (x,y,z) = C, 



(yz)%\ -y)(l- z)] a [(x + y + z){l- x)] a+p - l l 



2 



4" 



[(x + y)(l-x-y)]^(y + z) 



2;i 
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The exponents a and (3 are independent of n, quark helicity, and the transverse shape, 
while the auxiliary functions C n are dependent on all of these factors. If the C n are made 
from a complete set of functions of x a , non-singular in the limit when any number of 
momenta vanish, then any boundstate wavefunction should be expressible in the form 

We have also examined the effect of higher orders in n > 5. The forms ()57|) remain 
invariant, although the exponents a and f3 are renormlised in general because of parton 
self-energy and vacuum polarization. Wavefunctions ip n analogous to (|57|) become in- 
creasing complicated at higher n and we have not found a simple general expression; no 
new exponents beyond a and f3 occur however. Of special interest are the limits 

tjjnixt, x 2 , . . .) ~ (1 - Xl ) a -( n ~ 2)/2 as x x -> 1 (58) 
ipn{xi,x 2 , ■ • •) ~ as x\ -> (59) 

which follow from comparing interactions (z) and (ii). These govern the behaviour of 
hadronic structure functions at large and small Bjorken x, as discussed in Section 

The forms ([570 demonstrate that lightcone wavefunctions become more singular as 
more lightcone momenta vanish, even though this is a 'high-energy' limit. Intuition from 
the free lightcone kinetic energy (sum of all terms (i) and (vii) in Table ITT1 of Appendix IXj) 
would have suggested the complete opposite. Rather than simply vanishing at the edges 
of phase space, as indicated by the free theory, the interactions have forced relations 
between wavefunctions with different numbers of partons to cancel divergences in the 
boundstate equations. This was noted for continuum QCD for a single vanishing x a 
in Ref. jiJ], but for the transverse lattice we are able to generalise it to any number of 
vanishing momenta. Our analysis rules out a number of naive models for the lightcone 
wavefunctions, including products of individual parton wavefunctions and functions of 
the free kinetic energy. 



5. APPLICATIONS 

Based on the previous analysis, we can make definite statements about the behaviour 
of quark parton distributions in mesons at the edges of phase space. Some other recent 
studies of these distributions in various models have been made in Refs. [10|. We use the 



generalised quark distributions which in lightcone gauge are defined as 



1 f + °° H?' -—+ 

-H{x^Q 2 ) = -j=== J ^e^ MP in |*(-z-/2)7 + *(^2)|Pout> , (60) 



Q 



Pin 



Pout 

~ Pout; 



P^ 



2P 

i Pin + Pout; 



(61) 
(62) 



(63) 
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Bars generally indicate an average over in and out states. In Eq. ()60|). we define the 
boundstates at non-zero transverse momentum P by boosting the Fock states with the 
Poincare generators M + = (M +1 ,M +2 ): 

exp [-1M+P/P+] b{(x 1 P + )a{ 2 (x 2 P + )---a{ n _ i (x n ^ 1 P + )dl l (x n P + )\0) 
= exp [iP.c] 6j i (x 1 P + )4 2 (a; 2 P + ) • ■ ■ al n _ i (x n _ 1 P + )<,(x n P + )|0) , (64) 



c = (c\c 2 ) = y^x a x a , (65) 

a=l 

x p+ i = Xp + aX p+ i , p = 1 to n — 2 , (66) 

where A p+ i is a transverse unit vector in the direction of the p th link. The boundaries of 
interest are typically at large and small x and at the transition point x — £. 

A. x — > 1 

We find from Eqs. (f67IJ) and (JHEJ) that each ^ n contributes ~ (1 — x) 2a to H(x, £, Q 2 ) 
as x — > 1. This is interesting because the conventional wisdom, based upon the free 
kinetic energy, is that ip2 should dominate in this region — all but one parton has small 
momentum, so the free kinetic energy is minimized by the minimum number of partons. 
Although the exponent is correctly given by the lowest Fock wavefunction, which governs 
the distribution amplitude, on the transverse lattice all higher Fock states contribute to 
the structure function with the same power. 

B. x -> 

We find from Eqs. (jSQj) and that ip n contributes ~ (ln^)™ -3 , n > 2, to the 
traditional quark distribution function 7i(x, 0, 0) as x — > 0. Thus, the sum over all Fock 
sectors n rises faster than any power of lnx, i.e. at least like a power of x. This is 
consistent with the Regge behaviour of structure functions in the small-x region. 

c. x^i+ 

Unfortunately, we are not able to explore the region x < £ with the strict large- N c 
wavefunctions. That region is governed by their 1/N C corrections, since it depends upon 
quark pair production. However, we find from Eq. (|57J) that ip2 contributes zero, while ^3 
and ^4 contribute finite constants to 7i(x,C,,Q 2 ) asi->^ from above, provided £ > 0. 
This is presumably true for higher n also and, provided these constants fall sufficiently 
rapidly with n, suggests that 7Y(£,£,Q 2 ) is finite. 
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6. CONCLUSIONS 



By demanding finiteness of the boundstate energy, we have determined the form (|57J) 
of meson wavefunctions at small lightcone momentum in a coarse transverse lattice gauge 
theory at large N c . The forms depend on two exponents a and (3, independent of the 
other details of the wavefunction such as spin and number of partons, which encode 
the leading non-analytic behaviour. As more momenta vanish, wavefunctions become 
more singular, contrary to the naive intuition from free lightcone kinetic energy, which 
diverges in these limits. The analytic forms lead to predictions for the behaviour near 
their endpoints of quark distributions. In particular, we see Regge behaviour at small x 
and contributions from all Fock states at large x. 

The use of coarse lattice wavefunctions in transverse position space simplified the 
analysis, allowing us to separate the lightcone momentum problem from the transverse 
structure. The large- N c limit also aided us by restricting the diagrammatic rules to planar 
ones and suppressing quark pair production. We would need to reformulate the problem 
at finite iV c in order to study the small lightcone momentum behaviour in baryons and 
in the sea quark part of wavefunctions. There is also the question of higher orders of the 
colour-dieletric expansion; that is, higher powers of dynamical quark and link fields in 
the lightcone hamiltonian. We have investigated some of these possible terms, but found 
none that altered the form (|57j). although the values of a and (3 are affected. Finally, it 
would be interesting to investigate similar questions in continuum lightcone QCD. 
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APPENDIX A: INTERACTIONS 



(i) 



(ii) 



(iii) 



(iv) 



✓ 

z p ' 



(vii) 



X 



(X) 



x a* 



z p 



x a 



z p 



(xiii) 



y^ 



xh 



xh 



h' 



(v) 



y^ < ^ 



(viii) 



h' 
z p 



» (xi) 
yX ^ 

xa — > — x 
z p ' 



(vi) 



x h 



yX < \ 

/ 

/ 

z p '■ 



h' 



(ix) 



(xii) 



x 

y *v 

z p /' 



FIG. 5: Planar diagram representation of the elementary interaction vertices following at large 
N c . Time x + is increasing to the right. Solid lines are quarks, chain lines are link fields. 
x,y,z represent lightcone momentum fractions, h, h! helicities, X,a,p,r E {±1,±2} transverse 
orientation of links. Vertical barred lines are x + -instantaneous interactions. Spectator partons 
in the hadron are not drawn. 
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TABLE II: Matrix elements of the invariant mass operator 2P + P~ (for zero tranverse momen- 
tum) in Fock space; see Fig. [S] for corresponding diagrams of interaction vertices. Momentum 
conserving delta functions are omitted for clarity. P stands for principal part. Perp(A, p) = 1 
if |A| ^ \p | and otherwise. 



1 f.,2 | (kI+k 3 s )N c rx dy\ 
x\^f + 7T Jo y ) 





2^y{ {x 1 +y + l)^S^^ W c5hh'+( Kx l +y I ) M/«aV^c C\'h } 


(ti) 


-G 2 N C p / 1 \ 
2n r\{x- y y) 


(Hi) 


-G 2 N c (y~z) , 
4n(y+z) 2 ^yz°-*P 


(iv) 




(v) 


As (f ) with p — > —p 


(vi) 


m? 

b 

X 


(viz) 


—G 2 N C p ( y+z \ 

47r \Vz~y(z-y)' 2 ) 


(viii) 


-G 2 N c ~p ( (y+z)(2x+y-z) \ 
8,r \\fxzy{x+y-z){z-y) 2 J 


(ix) 


. , 1 T - r {2Ai(5 <tA SprS-pcj + \2(5<t\$pt5 P cj + 5-ap5-\ T 5- a \) 

47r-y/ xyz[x-\-z y) 


(x) 


+\4(5 a \5 pT Perp(a,p) + <5_ CTp 5_A r Perp((7, A)) - /3<5A /9 <5 CTT Perp(cr, p)} 




4^xyz(x + z+y) { 2X l S -^XrS Xp + A 2 <W£-a/|A| H 


(xi) 


+A 4 (<5_a ct 5 pr Perp(A, p) + (5_ CTp £vrPerp(<7, A)) - (36-\ p 6 aT Perp(a, A)} 




-C? 2 Af c (y~z)(2x+y+z) g 


(xii) 


—G 2 N C {x-z)(2y-x-z) g g 
87r \J xzy(x— y+z)(z+x) 2 °^ T 


(xiii) 
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